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1.  INTRODUCTION 

The  problem  of  determining  the  thermodynamics  of  a  Heisenberg 
ferromagnet,  as  a  function  of  the  temperature  and  the  magnetic 
field,  has  been  rigorously  studied  in  the  low  temperature  region^” 

m 

where  a  series  expansion  In  powers  of  (where  Tc  Is  the  Curie 

Temperature)  Is  valid.  Similarly,  at  high  temperatures  the  thermo 

dynamic  perturbation  theory^ “^hfti#  t>»ei»used  to  evaluate  the  free 

Tc  / 

energy  as  a  power  expansion  in  /^.  These  methods,  however,  are 

applicable  only  far  from  the  transition  region  and  therefore 

approximate  theories,  such  as  the  Weiss  molecular  fleld^’^theory 

and  the  various  cluster  theories, have  been  proposed  which 

give  qualilifctatlvely  good  results  near  and  beyond  the  transition 

temperature.  The  Curie  temperatures  can  also  be  estimated  from 

these  methods.  The  best  estimates,  however,  are  obtained  by  extra 
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polatlng  the  high  temperature  results  to  the  transition  region. 

All  these  methods  suffer  from  the  drawback  that  their  appli¬ 
cability  Is  limited  to  a  particular  range  of  temperatures. 

Recently  the  technique  of  double  tlrae-temperature  dependent 
Green  functions  '  has  been  successfully  applied  to  the  ferro¬ 
magnetic  problem.  The  virtue  of  this  method  Is  that  It  provides 
/  the  temperature  and  the  magnetic  field  dependence  of  the  magneti¬ 

zation  over  the  entire  temperature  range  with  reasonable  accuracy. 
Furthermore,  the  Curie  temperatxjres  provided  by  this  theory  seem 
to  be  In.  close  agreement  with  the  best  estimates  to  date  (compare 
references  8  and  21).  However  these  approximate  Green  function 
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theories  suffer  from  the  weakness  that  they  lack  a  detailed  agree¬ 
ment  with  the  exact  low  temperature  and  the  high  temperature 
expansions.  For  example,  the  low  temperature  result  for  the  raagnetl- 
zation  differs  from  the  spin  wave  theory  result  in  the  order  (T/Tcy. 
Similarly,  at  high  temperatures,  the  susceptibility  disagrees  with 
the  exact  result  In  the  order  (Tc/T)^. 

The  problem  of  Improving  on  the  Bogolyubov  and  Tyabllkov^"^”^®, 

21 

and  Tahlr-Khell  and  ter  Haar,  approximation  has  recently  been 

23 

studied  by  Callen  .  In  his  method  the  higher  order  Green  functions, 

-r 

occurlng  In  the  equations  of  motion  of  the  Green  functions  of  the 

spin  operators,  are  decoupled  by  a  plausible  physical  criterion. 

Callen' s  Ingenious  approximation  consists  essentially  in  taking 

into  account  the  fluctuations  of  the  operator  l.e.  the  z-com- 

g 

ponent  of  the  spin  operator  referred  to  lattice  site  g,  around 

z 

Its  statistical  average  <S>.  These  fluctuations  were  neglected 
in  the  earlier  work.  Callen' s  theory  successfully  predicts  the 
correct  spin  wave  energies  at  low  temperatures  and  also  leads  to 
an  accurate  estimate  of  the  Curie  temperatures  In  the  limit  of 
large  spin  values.  For  low  spins,  however,  the  results  are  less 
accurate  and  for  the  particular  case  of  S  =  1/2  the  expression 
for  the  low  temperature  magnetization  is  still  found  to  contain 

■5 

the  anomalous  (T/Tc)  term. 

At  high  temperatures,  Callen' s  theory  behaves  rather  similarly 
to  the  random  phase  theory*,  and  the  susceptibility  agrees  with  the 

18 

♦Here  and  henceforth  we  shall  refer  to  the  Tyabllkov’'*  ,  and  Tahir- 

21 

Kheli  and  ter  Haar,  approximation  as  the 'random  phase'  approxl- 

24 

matlon  (R.P.A.)  because  of  Its  equivalence  to  Engelert's  approximation. 
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2 

exact  results  up  to  order  (Tc/^)  . 

Prom  a  formal  point  of  view  both  the  random  phase  approxl- 
21  23 

mation  and  the  Callen  approximation  ^  assume  that  the  spin  Green 
function  obeys  the  Dyson  equation,  the  higher  order  Green  function 
being  written  as  a  product  of  a  mass  operator  and  a  lower  order  Green 
function.  These  theories  differ  in  the  form  that  they  postulate  for 
the  mass  operator. 

In  a  recent  study  of  the  problem,  Wortis  has  shown  that  the 

Green  functions  of  a  Heisenberg  ferromagnet,  in  common  with  those 

23 

pertaining  to  condensed  Bose  systems  do  not  obey  a  Dyson  equation 
with  a  simple  mass  operator;  alternatively.  If  an  "effective  mass 
operator"  Is  defined  by  a  Dyson  equation.  It  Is  found  to  possess  an 
anomalous  structure.  It  Is  our  purpose  here  to  investigate  the  form 
of  the  effective  mass  operator.  We  find  that  a  particularly  convenient 
representation  of  the  "effective  mass  operator"  is  such  that  the 
higher  order  Green  function  is  the  sum  of  an  anomalous  additive  term 
plus  the  product  of  a  simple  mass  operator  and  a  lower  order  Green 
function.  Knowing  the  detailed  form  of  the  spin  wave  dispersion  law, 

h_  27  2T 

’  ’  ' ’  and  keeping  In  mind  the  results  of  the  R. P.A.  and  the 

23 

Callen  theory  ,  we  Infer  the  form  of  the  anomalous  term  and  of  the 
remaining  mass  operator.  Similarly,  at  high  temperatures,  we  invoke 

the  series  expansions  for  the  susceptibility  and  the  magnetic  energy, 

8  9 

derived  by  Rushbrooke  and  Wood  and  Domb  and  Sykes  . 

We  find  that  the  mass  operator  Is  almost  exactly  that  proposed 

23 

by  Callen  At  low  temperatures,  the  anomalous  term  Is  quite  Important 
for  spin  ^  ,  but  it  becomes  Insignificant  for  higher  spins.  Again,  at 
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high  temperatures j  we  find  that  the  anomalous  term  contributes  only 

ij, 

weakly,  l.e.  In  the  order  (Tc/^)  ,  to  the  susceptibility  and  Is  less 

Important  the  larger  the  spin  S.  Near  the  Curie  point,  however,  the 

anomalous  term  makes  a  significant  contribution  even  for  moderately 

large  S,  l.e.  S;«^10.  For  S»l,  the  anomalous  term  Is  again  small 

23 

and  our  results  are  Identical  with  those  of  Callen 

The  Green  function  theory  here  obtained  provides  an  inter¬ 
polation  scheme  between  high  and  low  temperatures.  The  results  for  the 

C/vtcaJL  uo..tU, 

estimates  of  Domb  and  Sykes  and  Rushbrooke  and  Wood  to  about  1^  for  all 

spins.  The  critical  behavior  of  the  susceptibility,  as  T  approaches 

Tc  from  above,  and  of  the  magnetization,  as  T  approaches  Tc  from  below, 

is  Investigated.  It  Is  found  that  within  the  callen  and  the  random 

21 

phase  approximations  ,  the  susceptibility  Just  above  the  Curie  -feem- 
perature  has  a  f orm  X =  const. (1  -  Tc/^)"  whereas  the  magnetization 
Just  below  the  Curie  temperature  approaches  zero  as  (1  - 
The  present  theory,  however,  can  be  set  to  achieve  agreement  of  X.  with 
the  high  temperature  series  result  of  Domb  and  Sykes,  l.e  X  =  const. 
•(1-Tc/^)”^^^.  If  this  Is  done  we  find  that  Just  below  the  Curie  tem¬ 
perature  the  magnetization,  M(T),  obeys  a  relation  of  the  form  M(T) 

=  const.  (1-T/^^)^/^. 
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2.  THE  OREEN  PUNCTIOM 


The  mathematics  of  the  retarded  and  advanced-double  time- 

temperature  dependent  Green  functions  has  been  given  by  Zubarev^^« 

20 

and  by  Bonch-Bruevich  and  lyabllkov  to  whom  we  refer  for  details. 
We  outline  here  those  features  of  this  technique  which  are  relevant 
to  the  present  work. 

The  Green  function  of  operators  A  and  B>  «/i(t);  B(t'»,  Is 
defined  as  follows: 


4LACt)  ' 

«i»r. 

«Aet)  j 


.1 

— 

(2.1) 


where  A(t)  Is  the  Heisenberg  operator  referred  to  time  t,  l.e. 

ACO  =  A 

(2.2) 

and  where  H  Is  the  system  Hamiltonian,  2irfi  the  Planck  constant, 

A  the  Schroedlnger  (time  Independent)  operator,  square  brackets 
denote  a  commutator,  and  single  pointed  brackets  denote  an  ensemble 
average  _  /  -pM  \ 

<■■■> ' 

®<*)  Is  the  step  function 


f  1  ,  *  >o 

^  o  ,  X  <0 


(2.4) 


and  p  >  .A  j  ■  (kg:  Boltzmann's  constant;  T:  absolute  teaqperature). 
^B 
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Differentiating  the  Qreen  functions  given  In  (2.1)  with  respect 
to  t  we  get 


cLv  •  * 

(2.5) 


where  6 ( *  * ’ )  denotes  the  Dirac  delta  function  and  the  superscripts 
ret.  and  adv.  have  been  dropped  because  (2.5)  Is  the  same  for  either 
of  the  two  cases.  The  Green  function  «[A(t),H];  B(t')»  In 
general  Involves  Green  functions  of  higher  order  than  the  original 
«A(t);  B(t' )»j except,  of  co\u?se,  for  the  trivial  cases  of  non 
Interacting  systems  where  exact  solutions  can  be  obtained.  One 
has  therefore  to  linearize  the  eq\iatlon  (2.5)  by  a  suitable  de¬ 
coupling  approximation.  Once  (2.5)  has  been  solved  for 
«A(t);  B(t'  )»,  the  spectral  theorems^^"^®  may  be  Invoked  to  get 
the  time  correlation  functions. 


—  Lii^ 
€-»to 


(2.6) 


-  \ 


where  «AjB»^gj  denotes  the  energy  Fourier  transform  of 
«A(t)j  B(t’)».  Equations  (2.5)  and  (2.6)  are  the  only  two 
eq\iatlons  required  for  our  calculations  here. 
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3.  CORRELATION  FUNCTION  OP  THE  TRANSVERSE  COMPONBMTS  OP  SPINS 

We  shall  assiime  the  ferromagnetic  spin  system  to  be  described 
by  a  Heisenberg  type  Interaction  with  Isotropic  exchange  and  In  the 
presence  of  a  spatially  homogeneous  time  Independent  magnetic  field, 
B,  directed  along  the  positive  z-axls.  The  Hamiltonian,  H,  therefore 


Is 


H 


T 


where  E^  =  +  ^  Is  the  magnetic  moment  per  Ion, 

cartesian  components  of  the  spin  operator  for  the  site  f  and  I(f-m) 
the  exchange  Integral  between  Ions  at  sites  f  and  m.  As  us\ial,  we 
assvune  that  the  self-exchange  terms  vanish  l.e.  I(f-f)  ■  0,  We 
shall  consider  here  the  following  Green  functions: 


x,y,z 


(5.1) 

the 


where 


and 


,  c7ct')> 

=  s  * 


± 

S. 


"A, 


(5.2) 

(5.3) 

(3.4) 


qV,  =  [s/^')]  S,c^) 

where  n  Is  a  positive  Integer  or  zero. 

From  equations  (2.5)  and  3.1)  we  find  the  equation  of  motion 
/  \  A  ■  ♦  i**’’/  a 


where 


Q(a)*  4  (3.6) 
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In  order  to  solve  Eq.  (5.5)  for  the  Oreen  function  It  must 

first  be  linearized.  The  simplest  approximation  for  this  purpose 
Is  the  R.P.A.  ’  which  neglects  the  dynamical  correlations 

between  S^(t)  (or  S?(t)  )  and  the  remaining  operators  and  re- 

Z 

places  It  by  the  statistical  average  <S>.  In  order  to  take  into 
account  the  fluctuations  of  around  Its  average Callen^^ 
has  recently  suggested  a  decoupling  scheme  of  the  following  type: 


(5.7) 

Callen  chooses  a  on  the  grounds  that  since  the  random  phase  ap- 

the 

proxlmatlon,  represented  here  by  ^  choice  a  =  0,  affords  a  reason¬ 
able  first  approximation,  the  additional  part  -  <S^>  Introduced 

S 

here  must  be  self  consistently  small  at  all  temperatures.  This 

1 


requirement  can  be  satisfied  If  a 


at  low  temperatures  and 


If  a  decreases  at  least  as  fast  as  <S  >  at  high  temperatiu'es. 

z  2 

Callen' s  choice  of  or  =  <S  >/2S  Incorporates  both  these  features. 

In  order  to  look  for  an  Improvement  over  these  approximations, 
we  proceed  as  follows.  Rather  than  decoupllm^the^reen  functions 


«S^(t)  st(t);  cj"^(t')»  and  «S^(t)  S  ■'■(t)jl3^  Separately,  we 


-gv-,  , -  -g 

notice  that  In  equation  (5.5)  the  relevant  expression,  P'”*,  to 


be  decoupled  Is  a  function  of  the  difference  of  these  Oreen 
functions,  l.e. 


«  +  fc  ♦  /  u. 

<  * 


(J.8)a 
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Because  of  the  translational  and  the  time  Invariance  of  the 
Hamiltonian  and  the  translational  Invariance  of  the  lattice, 
must  be  a  function  of  the  vector  spatial  separation  (g-i)  and  the 


time  difference  (t-f),  l.e. 

f"'  h 

In  the  R.P.A.  the  function 


takes  the  following  simple  form' 
(*».>  Jus.) 


(5.8)b 

21 

1  . 

9) 


p'cvtjt-tO  <S*> 

In  view  of  the  fact  that  the  R.P.A.  constitutes  a  reasonable 
first  approximation,  we  propose  to  Investigate  a  decoupling  scheme 


of  the  form: 
I*'-) 


(5.10) 


where  represents  the  corrections  to  the  R.P.A. 

Introducing  (5.10)  Into  (5.5)  we  get  the  equation  of  motion 
in  the  following  form;  —  A 

The  translational  Invariance  dictates  Fourier  transformation 
with  respect  to  the  Inverse  lattice,  l.e. 


=  4- G  et 

3>1  ^ 


=  4:L 


M 


V  >«• 

e. 


(5.12) 

(5.15) 


(5.14) 


where  N  Is  the  total  number  of  sites  In  the  lattice  and  the  Inverse 
lattice  s\ms  are  restricted  to  the  first  Brlllouln  zone.  Fourier 


transformation  with  respect  to  the  energy  E  Is  also  possible, 
CA)  .  -l/RUi-'C;  ,, 

.<i= 


\  A  CK)Ot  .Je 


(5 

(5 


Using  (5.12)  -  (3.17)>  equation  (5*11)  takes  the  form 

C»V)  (»l) 

+  A  tKJE) 

where  Ej^  Is  the  elementary  excitation  energy  obtained  In  the  R 
(compare  equation  (3.11)  of  reference  21.),  l.e. 

Ek  =  3(0>k) 


Here  we  have  used  the  notation 

^  ICf)  ^ 

3Ck)  -  W)  - 


O 

(5 


It  Is  clear  that  In  a  simple  mass  operatorj[approxlmatlon,  the 
correction  will  be  of  the  form: 


.15) 

.16) 

.17) 

.18) 

.P.A. 

.19) 

.20) 

.21) 
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where  In  general^  would  be  complex  and  would  not  Involve 

terms  proportional  to  the  inverse  of  Q^”^(E).  Recently,  Wortls^ 
has  carried  out  a  careful  diagrammatic  study  of  the  structure  of 
the  Green  function  approximations  In  the  context  of  thez*modynamlc 
perturbation  theory  relevant  to  the  low  temperatures.  He  observes 
that  the  Green  function  pertaining  to  the  spin  problems  does  not. 

In  general,  have  the  structure  common  to  Green  functions  for 
normal  particle  systems  and  consequently,  the  mass  operator  -  the 
generalized  analogue  of  defined  In  (^.22)  -  has  an  ancMsalous 

structure  unique  to  the  spin  systems.  In  fact,  Wortls  notes  that 
whereas  In  the  case  of  the  normal  particle  systems  the  mass  operator 
Is  a  functional  of  the  Green  function  for  the  spin  systems  an 
effective  mass  operator  must  be  considered  to  be  a  functional  also 
of  the  Inverse  of  the  Green  function.  (Compare  argument  leading  to 
equation  (5-5)  of  reference  5).  For  the  present  purposes.  It  Is 
convenient  to  extract  the  part  proportional  to  the  Inverse  of 
G^*'^(E)  from  the  generalized,  effective  mass  operator,  l.e. 

Here  Is  a  function  of  the  system  temperature  and  of  the 

variables  n,  k  and  E,  and  M(k;E)  Is  the  analogue  of  the  us\ial  mass 
operator.  In  order  to  gain  some  Insight  Into  the  form  of  we 

Fourier  transform  (3.25)  as  follows: 


(5.24) 
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A  convenient  starting  approximation  for  can  now  be  arrived  at 
If  we  refer  to  equations  (3.10)  and ^.11).  It  Is  reasonable  to 
assume  that  because  of  the  discontinuous  nature  of  the  Green  function 
In  the  time  variable  (t-f)  -  refer  equations  (2.1)  -  (2,4)  - 
will  also  be  a  discontinuous  function  of  (t..t').  The  first  term  on 
the  r.h.s.  of  equation  (3.24)  already  Incorporates  this  discontinuity 
because  of  the  presence  of  the  Green  function  A  delta  function 

dependence  would  therefore  seem  to  suggest  Itself  for  the  remaining 
term  l.e. 

^  V  (3.25) 

where  y(T)  is  an  arbitrary  function  of  the  temperature  T,  R(g-.t)  Is 
a  suitable  function  of  the  vector  g-,^  and  the  temperature.  Here 
Q(n)  has  been  Introduced  for  convenience  In  later  calculations 
and  Is  the  same  as  In  (3.6).  Equations  (3.25),  (5.10)  and  (3.8)^^^ 
contain  the  essentials  of  the  form  of  the  basic  approximation  of 
our  theory.  Introducing  these  Into  the  equation  of  motion  (3.5) 
and  carrying  out  the  Fourier  transformations  described  In  (3.12)  - 
(3.17)  we  finally  get; 


(3.26) 


(compare,  equation  (3.l8). 

In  view  of  the  fact  that  the  Fourier  transform  of  the  Green 
function,  G^^^(E),  has  poles  at  B  -  +  I*^(E),  we  nay  recognize 
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it  to  be  the  analogue  of  the  elementary  excitations  of  the  system 
(see  ref.  20  for  a  discussion  of  this  point).  In  generak 
which  represents  the  shift  of  the  true  energy  spectrum  from  that 
obtained  In  the  R.P.A.^  will  be  complex.  For  the  present  p\irposes, 
however,  the  Imaginary  part  of  f^(E)  will  be  ignored.  It  should  be 
emphasized  that  strictly  speaking  the  damping  is  a  very  Important 
parameter  of  the  system  excitation.  It  determines.  In  fact,  the 
limits  of  the  applicability  of  the  concept  of  quasl-statlonarlty 
with  regard  to  the  elementary  excitation.  One  cannot,  therefore, 

/V 

claim  with  any  certainty  that  the  excitation  thus  chosen  re¬ 
presents  the  true  elementary  excitation  of  the  system  until  it 
can  be  ascertained  that  the  associated  damping  can  be  neglected. 

At  low  temperatures,  the  applicability  of  this  concept  Is  not  In 
doubt  because  the  associated  damping^  is  much  smaller  than 
Bearing  in  mind  this  conditional  interpretation  of  we  may 
proceed  as  follows: 

Prom  equations  (5.26)  and  (2.6)  and  the  Identity  that  for 
real  E  and  Ej^ 


r  ' 

1  , 


(5.27) 


we  get  the  following  expression  for  the  static  correlation  function 
of  and  Sg  :  ^ 


where 


K 

'♦v  _  4 


“  <C(\)  ^  ) 


(3-29) 
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$Ck) 


-  \) 


and 


(5.30) 


(3.51) 


Considering  the  fact  that 

*  5CS  +  I^^  - 


(3.32) 


when  £=  g,  equation  (3.31)  contains  both  on  the  left  and  right- 

^  z 

hand  sides  a  stun  of  averages  of  powers^S  .  We  can  therefore  write 

z 

down  2S  Independent,  simultaneous  linear  equations  In  <S>, 

<(S^)^, . .  <(S^)^^>  by  putting  n  in  equation  (3.28)  equal  to 

1,  2,  5,  2S  consecutively.  The  equations  with  n>2S  are 

not  Independent  of  the  earlier  ones  because  of  the  following 
operator  relation  satisfied  by  the  spin  operators: 

rT”7CS-r'l^')  =  o  (3.35) 

where  r  takes  on  Integral  or  half -odd-integral  values  according 
as  S  Is  Integral  or  half-odd-integral , 

The  results  for  the  average  <S^>  can  be  written  as  follows:* 

<1  > 1“’'! 

where 

5  S  §  +  ^  7J  ^  (3.35) 


20 

♦similar  expressions  were  empirically  arrived  at  after  solving 
explicitly  the  2S  simultaneous  equations  for  several  different 
values  of  S.  Callen^^  has  since  derived  these  results  by  a  much 
more  elegant  method. 
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The  expression  for  the  average  <(S  )  >  can  also  now  be  easily 
obtained  from  the  relation  (3.32). 

On  A*  proceeding  from  equations  (3.28)  and  (3.29)  to  equations 
(3.3^),  (5.35)  we  made  the  following  choice  for  the  function  R(lC); 

I  / 

X  $(|C) 


The  reasons  for  this  choice  become  clear  If  we  study  the  form  of 
the  transverse  correlation  function  of  the  spins,  L^°^(^-g). 
Putting  n»0  In  (3.2j)  we  get; 

CU-^)  =  <S^S.>  -  (3.37) 

At  low  temperatures,  where  It  has  already  been  arguedl  that  eT  of 
equation  (3.31)  will  represent  the  true  spin  wave  dispersion  law, 
the  form  of  the  correlation  function  -g)  Is  known  from  the 

spin  wave  theory.  It  Is  known  that  g,  this  correlation 

function  Is,  to  a  good  approximation,  given  by  the  first  term  on 
the  right-hand  side  of  (3.37).  Therefore,  a  convenient  choice 
for  R(k)  Is  given  In  (3.36)  because  then  the  remaining  term  In 
(3.37)  contributes  only  when  /  ■  g,  l.e.  L^°^(jt-g)  ^ 

-  +a,<§>-Y(T)  y  (5-58) 

Having  thus  specified,  up  to  an  arbitrary  temperature 
dependent  parameter  y(T),  the  form  of  the  function  of 

equation  (3.24)  we  Investigate  the  function  >^(^)*  Once  again 
we  Invoke  a  result  known  from  the  rigorous  spin  wave  theory  that 
at  low  temperatures 


fs$ 


E^-E, 


(J.39) 


(compare  equation  (3.7)  of  reference  4),  This  stiggeete  the 
following  form  for  ^^(E^^); 

•L 

^  T^O/jfO  (3.40) 

'  NJ  K' 

where  (r(T)  Is  to  be  a  temperature  dependent  parameter  such  that  In 
the  limit  of  low  temperatures  cr*(T)  approaches  unity.  Equations 
(3.28)  -(3.31),  (3.34)-  (3.35)  and  (3.4o)  constitute  a  set  of 
coupled  equations  which  must  be  solved  self  consistently  to  determine 
the  magnetization  and  the  transverse  correlation  function. 

The  formalism  developed  so  far  Is  valid  irrespective  of  the 
lattice  structure  and  the  spatial  dependence  of  the  exchange  Integral 
I(g-f).  In  order  to  facilitate  the  comparison  with  the  results  of 
other  theories  we  shall,  In  what  follows,  restrict  consideration 
to  lattices  of  cubic  symmetry  with  nearest  neighbor  exchange  Inter¬ 
action. 


~  ^  ^  ®  nearest  neighbors, 

=  (3.41) 

with  these  simplifying  assumptions,  the  elementary  excitation  energies 
of  equations  (3.31)  and  (3.40)^ take  the  following  simple  form; 

£.  t 


(3.42) 
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4.  THE  LOW  TEMPERATURE  REGION 


In  this  section  we  analyze  the  results  appropriate  to  tem^ 
peratures  low  compared  to  the  Curie  temperature. 

At  these  temperatures  the  magnetization  Is  close  to  the 
saturation  value  and  the  average  <i>  Is  little  different  from  S. 

The  function  $  therefore  Is  small  compared  with  unity  and  equation 
(5.54)  can  be  expanded  In  powers  of  %  to  give 


<S>  -(lS+i)  $  f  o($)  J 


(4.1) 


A  A  a,  V  A.*'  a.aS«-i  • 


AStA 


»St5. 


<(l) 


So  far  the  arbitrary  function  Y'(T)  remains  unspecified.  It  Is 
clear  from  the  definition  of  $  (see  eqns.  (5.50,  (5.42)  with 
(r(T)  =  1)  and  eqn.  (5.55)  that  at  low  temperatxires  It  Is  the 
strict  analogue  of  the  thermodynamic  average  of  the  number  of 
spin  waves  excited  per  lattice  site.  Olils  requires 


<!> 

and  therefore  y'(T)  Is  completely  specified: 

.AS 


(4.5) 


rex'!  =  V(T)A  =  fcstOi  —  -1 


As  expected,  y'(T)  Is  small  and  decreases  rapidly  with  the  Inoz^ase 
In  S.  In  fact,  y'(T)  Is  significant  only  for  S«^l.  Thus  the  slaqple 
me^ss  operator  approximation,  vdilch  Ignores  y'«  Is  sufficient  for 
large  S. 
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A/  .. 

The  calculation  of  ^  Involves  an  Integration  (or  rather  a 
summation)  over  the  Inverse  lattice  vector  k'.  The  Integrand, 

9(k')  J(k'),  on  the  other  hand,  depends  on  the  form  of  Ej^. 

Therefore  an  Iteration  process  has  to  be  used.  We  first  calculate 
<§>,  and  the  Integral  over  k'  In  equation  (5.^2),  In  the  R.P.A. 

The  second  Iteration  Is  obtained  by  Introducing  these  results  Into 
the  expression  for  Ej^  and  then 
At  low  ten^eratures,  this  Iteration  process  converges  very  fast  and 
It  t\irns  out  that  no  further  Iteration  Is  necessary  beyond  the  first 
Iteration  cycle  because  the  terms  not  Included  In  the  first  cycle 
contribute  In  a  higher  order  In  the  ratio  (T/Tc)  than  the  ones 
retained. 


a,  results  for  the  spin  wave  energies,  and  the  nagnetl- 
zatlon  M(T)  are  found  to  be  as  follows: 


T«^Tc  * 

(4.5) 

ncT) 

=  Mte)  ^ 

=•  Cf^«')/s)Ls  9^- -<a  9^*- 

(4.6) 

where 

(4.7) 

i/  ml,  for  slBq}le  cubic  lattice 

■  (5/4) (2)^'^,  body-centered  cubic 

■  (2)^/^,  face-centered  cubic 

(4.8) 

(4.9) 


(4.10) 


oCo  ’Z 

<^1  «  (3/i,)’TTV.25y,^ 

(kl.  =  TT^  t>S  )/^  Z 

®<^  =  Zs4  '  (!^s') 


and 


C>5 


simple  cubic 


=  — £  ^  body-centered  cubic 

288 

=  ^Vl6,  face-centered  cubic 


(4.11) 


These  results  are  the  same  as  obtained  by  Pyson^ (neglecting 
the  small  corrections  arising  from  the  2nd,  and  the  higher.  Born  ^ 
approximation  spin  wave  scattering).  The  results  for  the  average  <(?)> 
are  obtained  in  a  similar  fashion  and  we  get; 

/(if)  (“-la) 

.  a,  + 

.  ”c<,  i 

+  0  Cs'^)] 

For  the  case  S  =  1/2,  the  r.h.s.  of  equation  (4.12)  Is  exactly 
equal  to  4 . 

Combining  equations  (5.58),  (4.4)-(4.ll),  we  can  find  the 
correlation  function  of  the  transverse  coni>onent8  of  spins  appropriate 
to  the  low  tenQ)eratures, 


(4.13) 
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It  will  be  noticed  that  (4.1?)  Is  actually  an  Improvenent  on 
the  spin  wave  result  In  that  for  g  >  It  also  Includes  the  tern 
(last  one  on  the  r.h.s.  of  4.1?)  which  helps  fulfill  the  stringent 
requirements  of  the  spin  kinematics.  For  example >  for  the  simplest 


case  s  »  1/2,  we  have  the  exact  equality 


-  <<i> 


(4.14) 


This  equality  Is  satisfied  by  (4.1?)  whereas  the  usual  spin 
wave  result,  comprising  only  the  first  term  on  the  r.h.s.  of  (4.1?), 
does  not  satisfy  It. 


5.  THE  HIQH  TEMPERATURE  REGION 


In  this  section  we  consider  the  high  temperattire  expansions 
for  the  susceptibility  and  the  magnetic  specific  heat. 

In  the  presence  of  a  small  magnetic  field  B,  the  magneti¬ 
zation,  M(T)  y  Is  small  and  Is  proportional  to  B.  In  the  limit 
BsO,  the  energy  ^(see  equation  5.42)  Is  therefore  proportional 
to  <r(T)/B  as  long  as  cr(T)  does  not  go  to  zero  faster  than  B  ,  It 
Is  clear  from  the  results  of  references  l8,  21,  23  that  the  correct 

high  temperature  behavior  of  the  zero  field  susceptibility^ 

^  z 

quires  the  proportionality  of  E^^  to<S  >.  ISierefore,  It  Is  con¬ 
venient  to  put  ^ 

<rCT)  =  ^  ^CT)•  (5.1) 


where  A(T)  Is  a  function  of  the  teiiq;>erature  and  Is  to  be  non- 
Infinite,  l.e.  small  compared  with  (l/B).  The  energies  ^  thus  are 
(see  3.42)  of  the  following  form; 

=  E.  +  <s>  TO)  Ll  +  X(t)3  (5-2) 

where 


XCT 


^  ^  T(K)/j,o;  (5.5) 


and  L^°^(k')  Is  the  Fourier  transform  of  the  transverse  correlation 
function  given  In  (3.38).  As  the  self  exchange  Integral 

I(f-f)  has  been  taken  to  zero,  the  stun  ^  J(k:' )  vanishes  and 
therefore  (5.3)  Is  Independent  of  any  explicit  dependence  on  y(T). 
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In  this  temperature  range  the  quantity  f Is  large  conqpared  with 
unity.  Therefore  equation  (5.3^)  can  be  expanded  In  Inverse 

A/ 

powers  of  $  to  give  .|  'va"*  T 

The  zero  field  susceptibility, 1®  therefore  given  by  the 
following 

r u  .  y  c \  >  ■»  ~7  (e  u .  i\  i  I 

(5.6) 


oil  owing  ^  -1 


where  we  have  asstimed  that  a  reasonable  choice  of  y*  will  be  such 
that 


yCT)  « 


-1 


B 


(5.7) 


/\/ 


The  quantity  ($.B"^)  can  be  expanded  In  Inverse  powers  of  y 
where 


n'  P  T^o)C 

Introducing  the  quantities  t^  and  t^, 
to  ■  tanh  (PE^2) 

t^  .  taK)v[<S>  X«t))/HkY  7(0) J 


(5.8) 

(5.9) 

(5.10) 


and  using  equation  (3.2)  we  can  write  equation  (3.30)  as  follows: 

§CK)  *  -('/*')  (!4)[ '+  ^ 

=  -(va)  i-  )2 /t.^  J 
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Expanding  the  hyperbolic  tangent  tj^  and  Ignoring  terms  pro¬ 
portional  to  B  In  comparison  with  unity  we  get: 

$  =  ^  ^  (5.12) 

Fcti  [(i+ x(T';)<«>/4.'rt,J 


where 


I 

=  -ii 

IM 


(5.15) 


Before  discussing  further  the  evaluation  of  $  ,  we  recall  that  the 
R. P.A.  (represented  here  by  putting  X  *  0  and  y  »  0)  yielded  the 
correct  high  temperature  expansion  for  the  susceptibility,  up  to, 
and  Including,  terms  proportional  to  (  1/y  ),  Therefore  y  must 
Involve  terms  which  decrease  with  temperature  at  least  as  fast 
as  l/^"^. 


In  order  to  determine  the  correct  expressions  for  the  two 

parameters  A(T)  and  y'(T),  we  need  also  to  consider  the  correlation 

function  A  general  expression  for  L^°^(t-g)  la  not 

available  but  the  system  energy  and  the  magnetic  specific  heat, 

8  Q 

for  which  exact  expansions  are  known,  depend  sensitively  on 
L^°\£-g).  The  computation  of  the  magnetic  energy  requires  know¬ 
ledge  of  the  correlation  functions  of  the  transverse  as  well  as  the 
longitudinal  components  of  spins.  In  the  paramagnetic  region, 
spatial  Isotropy  requires  that  In  the  limit  B  *  0  the  longitudinal 
correlation  be  equal  to  the  transverse  one,  l.e. 


(5.14) 
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Thus  the  magnetic  energy,  <I>,  takes  the  following  simple  form 

/H>  =  (- Vz')Z  IC4-)<Ci>  =  (5.15). 

Combining  eqns.  (3.38)  and  (5.2)  we  get  the  average  magnetic  energy 
per  Ion 

C  _  (5.15)b 

^  (T>.Tcia-o') 

Thus  the  magnetic  energy  Is  explicitly  Independent  of  the  parameter 
Y'(T)  and  depends  only  the  parameter  A(T)  (or  X(T)  ). 

The  procedure  for  the  determination  of  A(T)  and  y*(T)  Is  now 
clear.  We  choose  A(T)  first  to  achieve  agreement  of  eqn.  (f*. 
with  the  result  of  the  known  high  temperature  expansion  for  the 
magnetic  energy.  y'(T)  Is  to  be  chosen  next  by  fitting  the  sus¬ 
ceptibility  expansion  of  eqns.  (5.6),  (5.12)  with  the  exact  one. 

Using  eqn.  (5.11),  eqn.  (5.15)h  can  be  expanded  In  Inverse 
powers  of  Y  bo  give 

e  =  f  c-o  [s(l  t  (5.16) 

where 

S  =  u 

f  ‘hm  O 

and  P(r)  Is  as  In  (5.13). 

The  stuns  F(r)  depend  only  on  the  crystal  structure  and  can 
be  calculated  easily  for  r>0.  (see  Appendix  A.)  The  quantity  S 
Is  proportional  to  the  zero  field  susceptibility  for  which  the 
exact  expansion  Is  known. 
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After  a  little  islgebra  one  oan  now  determine  the  high  tmi- 
perature  expansion  for  X(T)  suoh  that  agreement  between  eqn.  (5*15)b 
and  the  exact  expansion  for  the  magnetic  energy  Is  obtained.  Below 
we  give  the  results  for  the  two  leading  terms  In  the  expansions  for 
X(T)  and  the  related  quantity  A(T).  (Of.  eqns.  3.1  and  3.3) 

X(Y)  »  Xi//y  Y  (5  18) 


where 

=  C-'AV^ 


(5.19)0. 


Xt  =  (»4s  i**)  ®  O'/I  5  >'■) 

'  *  (5.»)k 


(iCp)*^-  29.630  f.c.c.j 


32.778  b.c.c;  69.136  simple  cubic 

(5.^S)^ 


and  where  z  denotes  the  number  of  nearest  neighbors.  Similar ly« 
we  have 

X  ^  ^,/y  ^  <5.20) 

A,  -  cvf V)  -  o/*c»»*))  ,5,21)^ 
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These  expansions  are  sufficient  to  achieve  agreement  for  the 

h 

magnetic  specific  heat  calculated  from  (3<15)  to  the  order  . 

The  susceptibility  expansion  (see  eqns.  5.6  and  5*12)  can  now 
be  similarly  derived  by  using  the  above  results  for  X(T).  An  In- 

e 

terestlng  fact  observed  Is  that  agreement  of  the  leading  three  terms 
In  the  l/y  expansion  for  Is  obtained  without  any  assistance  from 
the  parameter  y'(T).  The  fourth  term  Is  In  agreement  with  the  exact 
results  only  for  large  values  of  S. 

This,  once  again,  demonstrates  that  the  simple  mass  operator 
type  approximation,  Is  sufficient  to  describe  the  behavior 

of  spin  Green  functions  In  the  limit  of  large  spins.  For  general 
S  the  Y'(T)  necessary  to  achieve  agreement  for  the  susceptibility 

k 

up  to  and  Including  the  term  proportional  to  (1/^  )  #  followss 

-  '»*•  [VsVstol- •} 

As  long  as  the  exact  results  for  the  susceptibility  and  the  magnetic 
specific  heat  are  available  for  comparison,  the  fore-going  procedure 
can  be  used  to  determine  X(T)  and  y'(T),  and  thus  the  dynamical 
spin  correlation  functions,  to  any  order  In  l/rf. 
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6.  Temperatures  Close  to  To 

Unlike  the  case  of  the  very  low  or  the  very  high  temperatures, 

no  rigorous  expansions  for  thermodynamic  functions  exist  for  the 

range  of  Intermediate  temperatures  In  the  vicinity  of  the  Curie 

point.  It  Is  therefore  difficult  to  construct  a  formalism  which 

relies  on  the  knowledge  of  these  expansions  for  the  determination 

of  the  arbitrary  constants  #(T)  and  V(T). 

A  number  of  results  are,  however,  known  for  the  transition 

temperature  Itself.  For  Instance,  the  high  temperature  series 

8  Q 

extrapolation  techniques  provide  fairly  accurate  estimates  for 

the  Curie  temperature,  Tc,  as  a  function  of  the  exchange  Integral, 

I.  Similarly,  the  magnetic  energy  at  the  critical  point,  <I>^,  can 

c 

also  be  estimated  by  these  methods.  It  turns  out  that  these 
estimates  can  be  used  to  get  Information  about  the  functions  o(T) 
and  y(T)  In  this  temperature  range.  This  Investigation  forms  the 
contents  of  the  present  section. 

In  the  absence  of  applied  magnetic  field,  B  «  0,  the 

magnetization  M(T),  and  therefore  <S^>,  Is  small  and  ^  ^belng, 
under  these  conditions,  proportional  to  the  Inverse  of  <S^>^ls 
large.  Eqn.  (5.?^)  can  therefore  again  be  expanded  In  Inverse 
powers  of  ^  .  The  most  convenient  expansion.  It  turns  out.  Is  the 
*■  following  (compare  eqns.  5.4  and  5»5): 


<!>•<  - 

(6.1) 

where 

C, 

1  s  QascstO  - 

(6.2) 

c  •- 

r  -Mis'  -AH i  +•*‘4  s  »4IS  1 

(6.3) 

*■ 
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The  energies  ^  can  once  again  be  written  in  the  form  (5«2).  However, 

the  parameters  A(T)  and  X(t),  defined  in  eqxiatlons  (5»1)  and  (5.5)* 

will  In  general  have  different  temperature  dependence  from  that 

determined  In  the  previous  section.  The  appropriate  expansion  for 

equation  (5.50)  Is  now  the  following  (compare  (5*11)  ): 

|  =  jf<3)  (6.^) 

9  y  4r 

where  ^ 


Combining  equations  (6.1)-(6.5)  we  get: 

UscstO-l;  -  .  a.TcT)l 

^  »  +  X<T) 


(6.6) 


where  we  have  used  the  notation: 

PCTJ  -  <S>  YCT)  (6.7) 

The  fummatlon  pC-1)  (see  equation  5*1?)  is  well  known'*^  and  has 
the  following  values: 

P(-1)  =  1.54466,  f.c.c.;  1.59520,  b.c.c.;  l.  51658,  s.c.  (6.8) 

The  Curie  temperature,  T  ,  Is  obtained  by  requiring  that  (since 
z 

B  a  <S>^  0  as  T  approached  T^  from  below.  Thus  (6.6)  gives: 

.Y  _  r  .  .  (6.ql 


We  find,  empirically,  that  an  extremely  good  fit  of  the  results 
following  from  (6.9)  with  those  known  from  references  8  and  9  can 
be  obtained  (refer  Tables  I,  II,  III)  if  we  have: 
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(6.10) 


Eqn.  (6.10)  contains  two  constants,  X(Tc)  and  ^(To).  In  order  to 

specify  these  completely,  we  also  need  the  expression  for  the  critical 

value  of  the  magnetic  energj^  ^^T»Tc*  eqns,  (5.15)^^  ex- 

z  ^ 

pending  ^(k)  In  powers  of  <S>  (cf.,  eqn.  5.16),  and  proceeding  to 
^  z 

the  limit  <S>  =  o  we  get: 


(6.11) 


If  the  results  for 


f^/KsTc) 


were  available  for  all  values  of  S  and 
z,  equations  (6.10)  and  (6.11)  would  determlnd  X(Tg)  and  r(Tjj). 

This,  however.  Is  not  the  case. 

It  Is  convenient,  at  this  point,  to  establish  a  correspondence 
between  the  present  work  and  that  of  ref.  2?.  Callen  achieves  a 
mass  operator  type  solution  which.  In  our  notation.  Is  equivalent 
to  putting  (at  all  temperatures): 


^CT) 


('41)  5  '^(•0  0 


(6.12) 


Col\U.V 

At  low  temperatures,  our  results  were  obtained  by  choosing 
A(T)  »(l/2sj^f.  eqns.  5.4o  and  5.l)  and  Y 
Encouraged  by  this  correspondence  we  postulate 


=^5-,^  ['41-3  (6.13) 

However,  unlike  Callen^we  retain  the  anomalous  contribution  of  the 
mass  operator^  y(T).  Inserting  eqn.  (5*1)  into  (5*3)  and  putting 
T  «  Tc  we  now  get 
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(6.14) 


The  quantity^8^  ^  can  easily  be  calculated  from  eqns.  (6.11) 
and  (6.l4)  and  the  results  are  listed  in  Tables  IV,  V,  and  VI. 

We  notice  that  these  results  are  in  good  agreement  with  those 
available  from  reference  9  and  are  a  considerable  Improvement  on 
those  following  from  the  R.P.A. 

The  parameter  ^Tc)  Is  now  completely  determined  and  may 
easily  be  obtained  with  the  use  of  equation  (6.10)  and  X(To)  given 
In  Tables  IV-VI. 

We  shall  consider  next  the  behavior  of  the  parallel  suscepti¬ 
bility  at  temperatures  Just  beyond  the  Curie  point.  It  Is  con¬ 
venient  here  to  recast  equation  (6.1)  Into  the  following  form: 


f  =■  -  1  -  <S>  \_SC5+ij-i.1  fo-a^ 


3  <§> 

Let  us  put 

4  s  > 


(6.115) 


(6.16) 


where  )C  Is  proportional  to  the  zero  field  susceptibility  )C  ; 

^  /u* )  •  ;c 

(6.17) 

Introducing  (6.17)  (6*  tf),  expanding  ^  ,  and  proceeding  to 


x'  =  (s 


the  limit  Eq  -  0  we  get: 


—  (6.18) 
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Eqns.  (6.9)  and  (6.10),  which  define  the  Curie  temperatiire,  Tc, 
can  be  rewritten  as  follows: 


t;,. V-  '‘4>  -""I  -  .-5^,1 1  Z  “■«> 


Subtracting  (6.18)  from  (6.|^)  and  carrying  out  the  summations  over 
k  (see  Appendix  B)  we  find  that  Just  beyond  Tc  the  susceptibility 
X  has  the  following  form: 


(6.20) 


where 


^  =  (^) 


(6.21) 


If  we  follow  the  R.P.A.  and  the  Callen  assiunptlon  and  put  P(T)  «  o, 
or  alternatively,  if  we  make  the  approximation  that  In  the  vlnlnlty 
of  Tc,  r(T)  is  equal  to  r(Tc),  we  find  that  (6.20)  leads  to  the 
following 

(6.22) 


,31 


This  result  Is  similar  to  that  of  the  spherical  model-'*  and  was 

24 

earlier  obtained  for  the  spin  1/2  case  by  Bnglert  and  by  Kawasaki 
and  Morl^  .  Within  the  above  approximation  we  can  also  xindj  Just 
below  the  Curie  temperature  from  eqn.  (6.6). 


where 


f  =  F(-0  IC - 1 

L  i  +  xCTcl  J'-  F(-.)+6c,  3 


(6.24) 


( compare  ref .  21 ) . 


The  above  result  la  similar  In  form  to  that  of  the  molecular  field 
approximation  (M.P.A.) 


(6.25) 


where  Is  the  Curie  temperature  in  the  M.P.A. 

g 

Recently  Domb  and  Sykes^  have  reported,  after  a  painstaking 
examination  of  the  systematica  of  the  extrapolation  of  high  tem¬ 
perature  series,  that  Just  beyond  Tc  the  zero  field  susceptibility, 
X  t  obeys  a  relation  of  the  form 

-by. 

(6.26) 


rVs 


where  A  Is  a  constant  Independent  of  the  temperature.  It  Is  there¬ 
fore  clear  that  the  approximation  r(T)  *  I1[Tc),  which  led  to  the 
result  (6.22^),  is  unsatisfactory.  In  order  to  secure  agreement  of 
our  result  of  eqn.  (6. to)  and  the  above  result,  (6.^4),  the  function 
At),  In  the  Immediate  vicinity  of  Tc,  has  to  satisfy  the  following 
relation  & 

3SCsti)[  (6-2»7)^ 


4/ 


w 
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which  Is  equivalent  to  putting: 

p<T)  -  tii— (6-7)6 

3S  CS+OF 

where  C  Is  a  constant  Independent  of  the  temperature  and  n  >  4/3. 

Having  thus  Inferred  the  form  of  r(T)  close  to  Tc,  we  utilize 
It  In  the  determination  of  the  magnetization  Just  below  the  Curie 
temperature.  Inserting  eqn.  (6.27)^  Into  (6.6)  we  get: 

P  3't%)  U-t  (6.28) 

Using  eqns.  (3.35),  (3.38),  (5.3),  (6.4)  and  (6.13), 

z 

X(T)  can  easily  be  expanded  In  powers  of  <S>  and  we’ get: 

)(.<.T)Li+  K<T)3  =  X(T)t  J  (6.29) 

(compare  eqn.  6.l4). 
and  therefor*’  *■ 

XCT)«X(Tc)  ^4 


Y 

^oO-Yr.) 


(6.30) 


Thus  to  the  leading  power  In  the  difference  (1-1/,^,^)  we  may  replace 
X(T)  by  X(Tc)  In  (6.28).  Eqn.  (6.28)  now  easily  leads  to  the  result 

<s>=  k .  cp-F^t* +»<?-?* 

where  K  Is  to  be  a  real  quantity  and  Is  related  to  the  constant  C 
of  eqn.  (6.27)^  by  the  relation: 


=r— 

scs 


5tOpc3'^•)[[i■XCT«)3  [f 


ya 

fO  tael'll 


(6.32) 
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APPEMDIX  A 

The  Inverse  lattice  sums  In  equation  (5*13)  have  to  be  re¬ 
stricted  to  the  first  Brlllouln  zone.  It  Is  convenient  therefore 
to  transform  the  function  J(o,k)/J(o)  to  a  system  of  coordinates 
such  that  the  Integration  limits  are  simplified  (the  boundary  of  the 
Brlllouln  zone  In  the  cases  of  b.c.c.  and  f.c.c.  lattices  Is  stmeirtiat 
complicated  being  respeclvely  similar  In  structure  to  f.c.c.  and 
b.c.c.  lattices).  This  Is  done  by  choosing  the  coordinate  axes  In 
the  direction  of  the  basis-vectors  of  the  reciprocal  lattice  space, 
(for  the  basis  vectors  see  reference  33).  It  turns  out  that 


(A.l) 


where,  since  the  Jacobian  Is  unity  for  the  above  transformation,  we 
use  the  following  ||||escrlptlon  for  changing  sums  Into  Integrals: 


tiT 


(A.2) 


The  sums  F(r)  of  equation  (S^13)  are  now  elementary.  The  results 
for  the  first  several  oases  az*e  tabulated  below: 


I 


•  0  •  0  • 

1.000  000 
1.083  333 
1.222  222 
1.4l4  931 
1.668  403 
1.995  049 
2.412  423 

2.943  934 

3.620  184 


b.0.0. 

1.000  000 
1.125  000 
1.375  000 
1.802  734 
2.513  672 
3.696  533 
5.684  326 
9.066  343 
14.892  007 


8.0. 

1.000  000 
1.166  667 
1.500  000 
2.069  444 
3.013  889 
^.581  533 

7.209  619 
11.670  664 
19.338  445 
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APPBMDIX  B 

In  the  vicinity  of  Tc,  the  difference  between  X(T)  and  X(Tc)  Is  of 
the  order  <S>^  ^cee  Eqn.  (6,  Therefore,  subtracting  eqn.  (6.20) 

from  eqn.  (6.21)  we  get  for  temperatures  Just  above  the  Curie  teaqpera* 


ture 


SCSt  l)/-^]f  Cfc  -  ^)  t  SS  CSfO 

Ks 


Cit 


_ _  i.y  - 

*•  k  LV-t 


(B-l) 


where 


V  = 


*-X-'[h-<IT)3<  (B-2) 

The  Inverse  lattice  sununatlon  In  (B-l)  can  be  approximately 
performed  as  follows.  As  the  susceptibility  X.  Is  largely  Is 
small  compared  with  unity,  ^Bie  dominant  contribution  to  the  siuomatlon 
therefore  comes  from  small  values  of  k.  As  a  rough  approximation  we 
may  therefore  extend  the  Integration  limits  to  the  whole  of  the 
k-space  and  also  use  the  long  wave  length  approximation  for  J(o,k), 


l.e. 


J(0,k)  -  I(4^a^)  +  O^ka)*"  (B-3) 

where  the  nearest  neighbor  distance  has  been  taken  to  be^^i^z^(z 
Is  the  coordination  number).  In  this  manner  we  have: 

m 

(y  t  T(0A)3T(eiK;  ^  t  1k^ 

3  3  3 

where  v  Is  the  volume  per  Ion  (v  ■  a"^,  s.o.;  av2,  b.o.o.;  a  /4, 
f.c.c.). 

l^erefore,  eqns.  (B— >1)  and  (6.20)  give 


(B-4) 


TABLE  I  THE  PACE- CENTERED  CUBIC  LATTICE 
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TABLE  II 
KgT^'l  Rushbrooke  I^T^' 
'5F~J  Wood  l  lh^  , 


3Vc  > 

2J(0)S(S+1) 


BODY-CBMTERBD  CUBIC  LATTICE 
I  cur.  R.P.A. 


R.P.A. 


Vc 

Jg? 


Callen 


1/2 

2.60 

iL60 

jl.87 

ir 

1 

7.55 

7.66 

7.66 

3/2 

14.7 

l4.8i 

14.36 

2 

23.9 

24.05 

22.97 

5/2 

35.2 

35.39 

33.50 

5 

48.5 

48.82 

45.94 

s  0.752«  Rushbrooke  and  Wood 
=  0.785,  (Ours  and  Callen' s) 
•  0.718,  Random  phase  Approx. 


5.7 

9.1 

16.6 
26.2 
57.9 

51.7 


Rushbrooke 

and 

Wood 


TABLE  III  THE  SIMPLE  CUBIC  LATTICE 

l/Vc\  . . 


3V0 

2J(0)S(S+1) 


Ours 


(¥=) .  (&) 


Callen 


/2 

1.7 

1.75 

1.98 

2.7 

1 

5.25 

5.28 

5.28 

6.5 

1/2 

10.2 

10.27 

9.89 

11.8 

2 

16.65 

16.73 

15.83 

18.5 

)/2 

24.75 

24.65 

23.08 

26.8 

3 

33.9 

34.04 

31.65 

36.4 

■  0.716,  Rushbrooke  and  Wood 
-  0.734«  (Ours  and  Callen' s) 
*  0.659,  Random  phase  approx. 


f  ^  f 


TABLE  IV  THE  PACE~CENTBRBD  LATTICE 


r— 1 

V  Voj 

-  0.517,  for  all  S 

R.P.A. 

• 

s 

H 

X(Tc) 

1/2 

0.453 

0.439 

0.195 

1 

0.450 

0.449 

0.150 

3/2 

0.456 

- 

0.130 

2 

0.463 

- 

0.116 

5/2 

0.465 

- 

0.113 

3 

0.467 

- 

0.108 

m 

0.476 

0.474 

0.085 

TABLE  V  THE  BODY-CENTERED  CUBIC  LATTICE 


(  -^1  »  0.590,  for  all 

j  R.P.A. 

S. 

s 

X(Tc) 

1/2 

0.487 

0.211 

1 

0.508 

O.I6O 

3/2 

0.517 

0.l4l 

t 


TABLE  VI  THE  SIMPLE  CUBIC  LATTICE 


6VkbT^^h.p.a.  ■ 


s 

X(Tc) 

1/2 

0.626 

0.257 

1 

0.651 

0.190 

5/2 

0.664 

0.167 

2 

0.671 

0.155 

5/2 

0.674 

0.149 

5 

0.678 

0.142 

0.695 

0.114 

K  K  t 
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